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ABSTRACT 

Multi-epoch Chandra X-ray observations of nearby massive early-type galaxies open up the study of an im- 
portant regime of low-mass X-ray binary (LMXB) behavior — long term variability. In a companion paper, we 
report on the detection of 158 X-ray sources down to a detection/completeness limit of 0.6/1.4 x 10 37 ergss -1 
using five Chandra observations of NGC 4697, one of the nearest (1 1 .3 Mpc), optically luminous (Mb < -20), 
elliptical (E6) galaxy. In this paper, we report on the variability of LMXB candidates measured on timescales 
from seconds to years. At timescales of seconds to hours, we detect five sources with significant variability. 
Approximately 7% of sources show variability between any two observations, and 16 ±4% of sources do not 
have a constant luminosity over all five observations. Among variable sources, we identify eleven transient 
candidates, with which we estimate that if all LMXBs in NGC 4697 are long-term transients then they are on 
for ~ lOOyr and have a 7% duty cycle. These numbers are consistent with those found for brighter LMXBs 
in M87 and NGC 1399, which suggests that there does not appear to be a measurable difference between the 
outburst durations of long-term transient neutron star LMXBs and black hole LMXBs. We discuss in detail 
a transient supersoft source, whose properties are not easily explained by standard explanations for supersoft 
sources. 

Subject headings: binaries: close — galaxies: elliptical and lenticular, cD — galaxies: star clusters — globular 
clusters: general — X-rays: binaries — X-rays: galaxies 



1. INTRODUCTION 

With the angular resolution of the Chandra X-Ray Observa- 
tory low-mass X-ray binaries (LMXBs) in nearby early-type 
galaxies can be detected. For early-type galaxies with a rel- 
atively low X-ray-to-optical luminosity ratio, known as X- 
ray faint galaxies, the emission from LMXBs is the domi- 
nant source (e.g,. ISarazin et al . 2000, 2001, hereafter Papers 
I and II). In typical Chandra observations of such galaxies, 
tens to a few hundred LMXBs are resolved at bright X-ray 
luminosities, Lx > 5-10 x 10 37 ergss _I . By stacking multi- 
epoch Chandra observations of early-type galaxies, larger 
numbers of fainter LMXB candidates can be studied, and 
brighter LMXB candidates can be studied in greater detail. 
Such observations open up the study of an important regime 
of behavior — long term variability. Low mass X-ray bina- 
ries exhibit a wide range of variability, including transient 
outbursts, type I bursts, and milder fluctuations. Our under- 
standing of LMXB variability in the Galaxy drives models for 
accretion and LMXB evolution. However, there are several 
limitations in studying Galactic LMXBs: source distances are 
known for only a small subset, it is difficult to observe the 
whole Galaxy at once, and the size of the observed sample is 
very limited. Observations of nearby early-type galaxies can 
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overcome these limitations of Galactic studies. 

Examples of extragalactic LMXB candidates with extreme 
behaviors, such as luminous supersoft sources (SSs) and flar- 
ing sources, have also been examined with Chandra in greater 
detail. In our Galaxy and M31, SSs have very soft X-ray 
spectra (< 75 eV) and are generally believed to be accreting 
white dwarfs (WDs). Since the bolometric luminosities of 
extragalactic SSs can exceed the Eddington luminosity for a 
Chandrasekhar mass WD, an alternate hypothesis that these 
sources contain intermediate-mass (~ 10 2 - 10 3 M W ) accreting 
BHs has been put forth (Swart zet al.l2002T) . Due to the degra- 
dation of the soft response of the Advanced CCD Imaging 
Spectrometer (ACIS) on Chandra, the study of SSs has been 
limited. Some LMXBs exhibit another extreme behavior, 
namely relatively short timescale (seconds to minutes) flares. 
Recently, XMM-Newton detected type I X-ray bursts in M3 1 
dPietsch & Haberl 2005) and we reported a more ex treme ex- 
ample of flares in NGC 4697 (ISivakoff et al.ll2005l hereafter 
iPaper llHi . 

In this paper, we report on the variability of LMXB can- 
didates measured from multi-epoch Chandra X-ray observa- 
tions of NGC 4697, one of the nearest (11. 3 Mpc), optically 
luminous (Mb < -20), elliptical (E6) galaxy. In a compan- 
ion paper (Sivakoff et al 2008b, hereafter Paper IV) we report 
on the properties of LMXB candidates excluding variability, 
as well as the connection to globular clusters (GCs) revealed 
by a Hubble Space Telescope Advance Camera for Surveys 
(HST-ACS) observation near the galaxy center. 

In § |2] we recap key aspects of the observations and 
data reduction of NGC 4697 discussed in Paper IV. Our ob- 
servations of NGC 4697 provide five approximately equal 
chances to detect variability within a given ~ 40 ks observa- 
tion (§|3]l, as well as variability between different observations 
on timescales of 1 1 d to 4.6 yr (§|4]l. We summarize our con- 
clusions in §|5] Unless otherwise noted, all errors refer to lcr 
confidence intervals, count rates are in the 0.3-6 keV band, 
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and fluxes and luminosities are in the 0.3-10keV band. 

2. CHANDRA OBSERVATIONS AND DATA REDUCTION 

Since detailed descriptions of the Chandra observation and 
data reduction are reported in Paper IV, we only recap key 
aspects here. This analysis covers five Chandra ACIS ob- 
servations (Observations 0784, 4727, 4728, 4729, and 4730) 
taken on 2000 January 15, 2003 December 26, 2004 January 
06, February 02, and August 18. The final flare-filtered live 
exposure times for the five observations were 37174, 39919, 
35601, 32038, and 40044s. There are differences in the detec- 
tor setup and instrument properties, such that the response of 
the Cycle- 1 observation (0784) differs most from the Cycle-5 
observations (4727, 4728, 4729, 4730). Key differences in- 
clude the gain files used (acisD1999-09-16gainN0005.fits for 
Observation 0784 versus acisD2000-01-29gain_ctiN0001.fits 
with time-dependent gain and charge transfer inefficiency cor- 
rections for the Cycle-5 observations), smaller background 
due to Very-Faint mode cleaning in Cycle-5 observations, and 
increasing quantum efficiency (QE) degradation with time. 

All Chandra observations were analyzed using CIAO 3.1 7 
with CALDB 2.28 and NASA's FTOOLS 5.3 8 . Source posi- 
tions and extraction regions were refined using ACIS Extract 
3.34 9 . All spectra were fit using XSPEcP. 

The finding charts for the sources are displayed in Figures 
1 and 2 of Paper IV. In that Paper, the source properties from 
stacking the multi-epoch observations are listed in Tables 1 
and 7, and the instantaneous luminosities for each observation 
are given in Table 6. 

Throughout Paper IV and this paper, we refer to the Anal- 
ysis Sample, which includes 126 sources that have pho- 
tometric count rates determined at the > 3a level when 
stacking all five observations. This corresponds to having 
at least 18 net counts, a minimum detected count rate of 
1.0 x 10~ 4 counts s -1 , and a minimum luminosity of 1.4 x 
10 37 ergss _1 . 

The spectral properties of sources can be crude ly char- 
acterized by hardness ratios or X-ray colors (e.g., iPaper It 
iPaper It) . We defined hardness ratios of H 2 \ = (M-S)/{M+ 
S), H 3 i = (H-S)/(H+S), and // 32 = (H-M)/(H+M), where 
S, M, and H are the total counts in the soft (0.3-1 keV), 
mediu m (1-2 keV), and hard (2-6 keV) bands (Sivakoff et al. 
l200l) . 

3. INTRAOBSERVATION VARIABILITY 

After applying barycentric corrections to the times in event 
files, we applied three techniques to searc h for intraobserva - 
tion variability: th e Raylei gh statistic (e.g jMuno et al.l2003l) . 
the K-S test (e.g., Paper nh, and a newly developed flare de- 
tection technique dPaper II J) . These three techniques consider 
all photons in the source apertures, and thus are likely to in- 
clude some background photons. 

We used the Rayleigh statistic to search for periodic sig- 
nals with frequencies between 5.0 x 10~ 5 Hz and 1 .0 x 10 _1 Hz 
(periods of 10-20, 000s), testing every 2.5 x 10" 5 Hz. TableQ] 
summarizes the results for > 2a level variable sources (i.e., 
the chance the observed source variability is due to a statisti- 
cal fluctuation is < 4.6%), listing the observation, probability 
that the detected periodicity is due to a statistical fluctuation, 
Pr, and detected period, tr. Between zero and two sources 



TABLE 1 

Possible Periodically Variable Sources 









Period tr 


Source 


Observation 


Probability Constant 


(sec) 


(1) 


(2) 


(3) 


(4) 


1 


0784 


1.1E-02 


21.93 


101 


4727 


6.7E-03 


14.98 


103 


4728 


4.4E-02 


20.02 


143 


0784 


4.6E-08 


1000.00 


155 


4728 


8.1E-03 


16.99 



NOTE. — The value of Pr indicates the probability that the 
observed periodicity is due to a statistical fluctuation. All periods, 
tr, are measured in seconds. No periodic variable sources were 
detected in Observations 4729 and 4730. 
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7 See http://asc. harvard. edu/ciao/ 

8 See http://heasarc.gsfc.nasa. gov/docs/sof tware/lheas 
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FIG. 1. — Binned lightcurve using 1350 s bins of Source 156 (the fore- 
ground star BD-05 3573) in Observation 4727. Both the K-S test and our 
flare detection technique classify this source as variable within this observa- 
tion. This star flared shortly after the observation began. The average count 
rate during the observations is indicated with the dotted line. The count rate 
during the flare determined by our flare detection algorithm is indicated by 
the dashed line and the extent of the flare is indicated by the gray region. 

appear periodic at the > 2a level in any given observation. 
Since the expected number of false detected sources at this 
level is ~ 7, we concentrate our discussion on the one pe- 
riodic source detected at the > 3a level (probability that the 
variability is due to a fluctuation is < 0.27%). The probability 
that the apparent periodicity among the 26 counts of Source 
143 in Observation 0784 is due to a statistical fluctuation is 
only 4.6 x 10~ 8 . Unfortunately, this source is located near an 
edge of the S3 chip, and the period of ~ 1000s is consistent 
with a periodicity induced by the yaw of the satellite aspect 
motion. Thus, it is likely that this periodicity is instrumen- 
tal, and we have no cases of detected periodic variations for 
sources in NGC 4697. 

By comparing the cumulative fraction of events received 
from a source to a constant rate, the K-S test can identify 
intraobservation variability. This method ha s been a pplied 
to LMXB candidates in previous papers (e.g.. lPaper It) . We 
summarize the results in Table [2] for sources which are vari- 
able at the > 2a level. We list the probability the detected 
variability is due to a statistical fluctuation for each obser- 
vation, Pks- Between one and nine sources appear variable 
in any given observation. Again, we concentrate on the one 
variable source detected at the > 3a level. Source 156 is 
the foreground star BD-05 3573 which undergoes a clear 
flare (Pks = 6.9 x 10~ 7 ) soon after Observation 4727 begins. 
We display the binned lightcurve of this source in Figure [JJ 
of We note that source 137 appears as a 2cr variable source in 
Observations 0784 and 4727. Following equation 3 of 



VARIABILITY IN THE LMXBS OF NGC 4697 



3 



TABLE 2 
KS Test for 
intraobservation 
Variability 



Source 
(1) 


Observation 
(2) 


(3) 


2 


0784 


2.4E-02 


10 


4727 


4.6E-02 


11 


0784 


1.4E-02 


12 


4728 


1.2E-02 


40 


0784 


2.9E-02 


42 


4730 


3.0E-02 


54 


4729 


1.2E-02 


65 


4727 


1.6E-02 


71 


4728 


5.0E-03 


75 


4727 


2.0E-02 


76 


4729 


4.3E-02 


78 


0784 


1.5E-02 


79 


0784 


1.5E-02 


83 


0784 


2.6E-02 


96 


0784 


3.5E-02 


100 


4729 


1.6E-02 


115 


4728 


3.9E-02 


130 


0784 


2.0E-02 


137 


0784 


3.7E-02 




4727 


2.0E-02 


156 


4727 


6.9E-07 


NOTE. — The value of P KS in- 
dicates the probability that the ob- 



served variability is due to a statis- 
tical fluctuation. 

iPaper III the joint probability that this source is constant is 
4.6 x 10~ 3 , which is still less than 3a significant. 

Finally, we applied a new method developed to search for 
flaring sources. This technique is based on the arrival times 
of individual events compared to the Poisson distrib ution at 
a constant rate, and was presented previously in IPaper Ilj . 
This flare detection uses the ONTIME values (37651, 40447, 
36072, 32462, and 40574 s, respectively for observations 
0784 and 4727^-730) as opposed to the live exposure times. 
IPaper IlJ gave the results for the most significantly flaring 
sources; note that Sources A, B, and C in that Paper cor- 
respond to Sources 71, 75, and 57 here, respectively. We 
summarize the results of these tests in Table [3] for sources 
which are variable at the > 2a level in at least one observa- 
tion. For each observation of a flare-candidate source, we list 
the observation number, the total number of photons in the 
observation, N, the number of photons in the flare, n, the du- 
ration of the flare, 8t, and the Modified Julian Date (MJD) 
of the arrival time of the first photon in the flare, to- The 
value of ^constant indicates the mathematical probability of the 
given flare occurring due to a statistical fluctuation in a con- 
stant rate source; however, we must account for the number 
of non-independent flare searches undertaken. After applying 
corrections to P con stant based on the number of photons in the 
observation, P' mn;i teni i s tne probability a flare in a given obser- 
vation is not real. The number of sources with a statistically 
significant flare ranges from one to seven in the different ob- 
servations. 

Given the number of sources expected to show a flare from 
a statistical fluctuation, we only consider the flaring sources 
in individual observations detected at the > 3a level, which 
are Sources 57, 155, and 156. These sources have flares rang- 
ing from a few photons over tens of seconds to nearly a hun- 
dred photons ov er thousan ds of seconds. Source 57 was dis- 
cussed earlier in lPaper ill] as Source C; note that the luminos- 



Src. 57 Obs. 0784 





i 














No 


a 






1 oooc 


Exposur 
Flare C 


e I 


c 


ne (s) 
idate 


30000 



21800 22200 



FIG. 2. — Impulse diagram indicating the time of arrival of photons in 
Source 57 (Observation 0784). This LM XB candidate undergoes short, bright 
flares in multiple observations I Paper III), which were detected only with our 
new flare detection technique. 
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FIG. 3. — Cumulative lightcurve of Source 155 (Observation 4727) over- 
laid by the expected distribution for a constant source (dotted line). Our flare 
detection technique classifies this source as variable within this observation; 
however, the K-S test misses the variability of this source. The flare inter- 
val, which lies between the dashed lines, has a clearly increased rate (larger 
slope). 

ity conve rsion used in this Paper is about 10% smaller than in 
IPaper H I due to updated calibrations. We display an impulse 
diagram of the photon arrival times for Observation 0784 of 
Source 57 in Figure [2] The observed count rate of the flare is 
over 100 times the average count rate. Source 57 has no opti- 
cal counterpart in the HST-ACS image, suggesting it is most 
likely a flaring Field-LMXB. 

A flare of several hundred seconds with 17 counts is 
detected in Source 155, but was missed by the K-S test. 
This source was just above the criteria used in IPaper IIll 
(^constantjoint < 0.32%) to identify a flaring source. The ob- 
served count rate of the flare is about five times the average. 
As shown in Figure [3] the flare is marked by a clearly in- 
creased rate (larger slope in the cumulative light curve) that 
straddles the cumulative light curve for a constant rate source. 
Since the K-S test looks for the largest deviation from a con- 
stant rate, it is not surprising that the K-S test fails to de- 
tect this type of flare. Follow-up of the optical counterpart 
to Source 155 is necessary to determine if this source is an 
flaring LMXB in a GC or a source unrelated to NGC 4697. 
As discussed earlier, Source 156 is the foreground star BD- 
05 3573 that undergoes a clear flare with 99 counts. In Fig- 
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TABLE 3 

Possible Flaring Sources in NGC 4697 



Source Observation N n At t ^constant Constant P constant,joint 

(s) (MJD) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



25 


4730 


21 


9 


2952 


53235.81643 


2.5E-03 


6.5E-03 


3.2E-02 


42 


4728 


26 


3 


25 


53010.59465 


3.9E-03 


1.0E-02 


5.2E-02 


56 


4727 


10 


5 


1919 


coaaa onoTi 

52999.89337 


8.7E-03 


2.1E-02 


1.1E-01 


57 


0784 


20 


4 


68 


51558.95770 


1.3E-04 


3.4E-04 






4727 


20 


3 


50 


COAAA A^in 

53000.06223 


5.5E-03 


1 .4E-02 






0784, 4727 














4.8E-05 


71 


0784 


14 


5 


1047 


51558.79423 


7.0E-03 


1.8E-02 






4727 


16 


4 


628 


52999.75882 


2.7E-02 


7.0E-02 






4728 


14 


4 


509 


coaia ^"7 a o a c 

53010.70295 


1.2E-02 


3.1E-02 






0784, 4727, 4728 














3.9E-04 


75 


4727 


10 


5 


1329 


CI AAA A A O O i~ 

53000.04326 


O 'TP AO 

2.2E-03 


5.5E-03 






4728 


9 


4 


1420 


C O A 1 A A A A 

53010.72444 


3.4E-02 


o nr ao 

8.0L-02 






4729 


6 


5 


5654 


c ~> r\ rani c 

53047.59915 


4.3E-02 


A AT 1 AO 

9.0E-02 






4727, 4728 














4.4E-03 




(Til A 

4727^4-729 














4.0E-04 


78 


0784 


44 


4 


52 


ri fro n^l A A 

51558.84209 


1.5E-03 


4. 1E-03 


O AT 1 AO 

2.0E-02 


83 


0784 


10 


8 


6910 


51558.83658 


8.5E-03 


2.1E-02 


1.1E-01 


90 


4727 


13 


3 


1 16 


52999.82502 


7.5E-03 


1.9E-02 


9.5E-02 


107 


0784 


2 


2 


480 


51559.09218 


2.5E-02 


5.2E-02 






4728 


4 


2 


63 


C O A 1 A z' A A O A 

53010.69939 


2. 1E-02 


4.2E-02 






0784, 4728 














2.2E-02 


115 


4728 


22 


8 


1857 


53010.37696 


2.7E-03 


6.9E-03 


3.5E-02 


118 


4728 


35 


11 


2667 


53010.42386 


9.0E-03 


2.5E-02 


1.2E-02 


124 


4728 


20 


4 


270 


53010.38518 


8.5E-03 


2.2E-02 


1.1E-01 


135 


0784 


4 


2 


78 


51558.70911 


2.4E-02 


5.0E-02 






4730 


6 


2 


57 


53235.77593 


4.1E-02 


8.5E-02 






0784, 4730 














4.3E-02 


137 


4727 


27 


12 


5529 


53000.01065 


2.4E-02 


6.5E-02 






4729 


18 


6 


1034 


53047.40640 


4.2E-03 


1.1E-02 






4727, 4729 














6.9E-03 


146 


0784 


13 


4 


697 


51559.01802 


1.9E-02 


4.9E-02 






4727 


4 


2 


22 


52999.76129 


6.5E-03 


1.3E-02 






0784, 4727 














6.6E-03 


155 


4728 


219 


17 


584 


53010.48733 


2.2E-04 


7.5E-04 


3.7E-03 


156 


4727 


289 


99 


7423 


52999.64230 


4.0E-06 


1.3E-05 


6.7E-05 



NOTE. — This table follows the nomenclature of Pacer III: N is the total number of photons during an 
observation. A flare candidate with n photons of observed duration St begins at to, where we list the time 
measured prior to barycentric corrections. The value of P CO nstant indicates the mathematical probability of 
the given flare occurring due to a statistical fluctuation; Constant ' s me same probability after accounting for 
the number of non-independent flare searches the technique examines. If the same type of flare occurs in 
multiple observations, one can calculate the joint probability the flare is not real, P C onstant,joint- 



ure[TJ we overlay the flare count rate determined by the flare 
detection algorithm for Source 156, which is about twice the 
average rate. 

If we consider searching for flares over multiple observa- 
tions, one can calculate the joint probability the flares are not 
real, P C onstant.joint- Thirteen sources in the Analysis Sample 
have a ^constantjoim that indicate a > 2a detection of flaring 
across multiple observations; however, the detected flares in 
two of those sources (137 and 146) have disparate enough 
durations that their separate flare detections probably should 
not be combined. Only 4 sources (57, 71, 75, and 156) have 
flaring which is significant at the > 3a level acro ss multipl e 
observations; these are the sources discussed in (Pa per II j) . 
plus the foreground star BD-05 3573. After correcting for 
the number of sources expected from statistical fluctuations, 
we find that 4.2^ '\% of the Analysis Sample sources exhibit 
flaring which is detectable at the > 2a level with our new al- 
gorithm. This is consistent with the 3.1*j'j% of the Analysis 
Sample sources that exhibit flaring which is detectable at the 
> 3(7 level. 



4. INTEROBSERVATION VARIABILITY 



Our observations of NGC 4697 are among the first designed 
to measure interobservation variability with strong coverage 
over a wide time baseline. Some variab ility results have been 
presented before (e.g., Centau rus A: iKraft et all 1200 ll M87 
Jorda n et alJ|2004Ulrwinf 2006. NGC 13 99lLoewenstein et al ' 
12005b llrwrnl I2006L and NGC 4636: iPosson-Brown et al 



2008); however, these results did not have the same qual- 
ity temporal coverage as our NGC 4697 observations. More 
recently, variability resu lts were discussed for NGC 3379 
dBrassington et alj 120081) . with a similar coverage of time 
baselines as in our NGC 4697 observations. We searched for 
luminosity variability between any two observations, as well 
as luminosity or spectral variability over the entire set of ob- 
servations. 

4.1. Luminosity Variability Between Two Observations 

We first searched for variability between pairs of observa- 
tions by calculating the value of \ 2 assuming a constant flux. 
Sources were considered possibly variable if \ 2 > 4, i.e., a 
2a detection of variability. This method differs slightly from 
IPosson-Brown et all d2008l) . which compares the individual 
luminosities to each other. The probability that the variability 
is due to a statistical fluctuation is listed in columns (2)-(ll) 
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FIG. 4. — Percentage of Analysis Sample sources that appeal' variable (at 
the > 2ct level) between two observations against the time period between the 
observations. Note that the error bars are not independent as the pairs contain 
the same five observations. The dotted horizontal line represents the percent- 
age of sources that would appear variable just due to statistical fluctuations. 
Although there is a rough tendency for more variability on longer timescales, 
a constant variability percentage of 11.3% is statistically consistent with the 
data. 

of Table H] for the ten pairs of observations. In Figure HI we 
display the percentage of variable sources as a function of the 
timescale between two observations for all of the Analysis 
Sample sources that are in the FOV for both observations in 
each pair. We also display the percentage of variable sources 
expected due to statistical fluctuations at the 2a level. There 
is a slight tendency for a higher percentage of variable sources 
at longer timescales; however, this is not statistically signifi- 
cant, based on a \ 2 test- (There are two caveats with this: the 
error bars in Figure [4] are not independent, and there may be 
a slight excess variability on the four longest timescales due 
to the improper correction of QE degradation for SSs. This 
latter effect comes about because we have assumed an aver- 
age spectrum when correcting for QE degradation; however 
SSs have a much softer spectrum that is more strongly af- 
fected by the QE degradation and the differences due to QE 
degradation will be most evident in comparisons to the first 
observation. We estimate that this latter effect is well within 
the displayed errors that arise from counting statistics.) If we 
treat the error bars as independent, we find that 11.3 ± 1.0% 
of the Analysis Sample sources are variable on average at 
the 2a level. Given the expected number of falsely identified 
variable sources, 6.7 ± 1.0% of the Analysis Sample sources 
are intrinsically variable between any pair of observations. If 
we only consider sources variable at the 3a level, the frac- 
tion of variable sources is 6.0 ± 0.7%. Thus, over timescales 
of ~ 10-2000 days, ~ 6% of sources are significantly vari- 
able. This number is roughly consistent with that found for 
Centaurus A (9. 8 +^%) for a variability timescale of 164d 
dKraft et al. 2001); although they derive a much higher per- 
centage when they include a generous definition of transient 
candidates (14.2+2 3%)- O ul derived variability percentage is 
lower than that found for a 1 1 1 5 d timescale in NGC 4636 
(23.8 + 3 Q%: lPosson-Brown et alj|2008h ; however, it is consis- 
tent with the variability for a day timescale (4.6+2 

4.2. Long-Term Luminosity Variability 

We can apply a similar x 2 technique to search for sources 
that show any luminosity variability over the entire period of 
five observations. We refer to these sources as Long-Term 
Luminosity Variables (LTLVs), and the probability Pl.au that 
the variability arises from statistical fluctuations is given in 
column 12 of Table [4] The difference in % 2 compared to a 



single merged luminosity and combined with the appropriate 
degrees of freedom (one less than the number of observations 
a source is in the FOV) is used to calculate the probability an 
observed LTLV is due to a statistical fluctuation. Sources 148, 
154, and 157 are excluded from these calculations as they are 
only in the FOV of one observation. From the Analysis Sam- 
ple, 26/124 sources (21.C£f;jj%) are LTLVs at the 2er signifi- 
cance level, while 11.3+1 9% of Analysis Sample sources are 
LTLVs at the > 3er significance level. Statistically correct- 
ing for falsely identified LTLVs, we find 16.4+j \% of the sig- 
nificantly detected sources show LTLV behavior at the > 2a 
significance level. Three of these sources are not associated 
with NGC 4697: Source 45 has an optical ID that is not a GC 
and thus is probably due to a background AGN, Source 1 17 is 
a known background AGN, and Source 156 is a known fore- 
ground star. Two sources (25 and 155) are associated with op- 
tical counterparts that may or may not be GCs in NGC 4697. 
Sources 39, 41, 51, 77, and 79 are GC-LMXBs that exhibit 
LTLV behavior. (Sources 98 and 114 are associated with 
Kavelaars GC candidates; new HST-ACS data will be used to 
test whether these sources are really GC-LMXBs.) After cor- 
recting for the likely percentage of falsely identified LTLVs, 
10.1+|;^% of GC-LMXBs and < 8.5% of Field-LMXBs ex- 
hibit LTLV behavior at the > 2a significance level. 

The recent analysis of sources in NGC 3379 found 
that 42 ±5% of their LMXBs exhibited LTLV behavior 
dBrassington et al.l [2008) using a % 2 test similar to the one 
used in this paper; however, they only required a source have 
a reduced x 2 > 1 -2. For sources in 2-5 observations, this cor- 
responds to a ~ 30% false detection rate. After correcting 
for these false detections, the rate for LTLVs in NGC 4697 
(~ 11%) and NGC 3379 (~ 12%) are comparable. 

Since some Galactic LMXBs are known to undergo tran- 
sitions between distinct luminosity states, we developed a 
method that groups observations together into possible lumi- 
nosity states. For every LTLV, we first considered all possible 
combinations of observation groupings that divide into two 
separate states. That is, we assumed that the luminosity had 
two distinct values, and that one or the other of the two values 
applied to each observation of that source. The two lumi- 
nosity values were allowed to vary, and each observation was 
assigned to one of the two luminosity states so as to minimize 
X 2 - Some sources still had a much larger x 2 than the number 
of dofs. Given the new x 2 , we could calculate the probability 
that any deviations from our two-state hypothesis were due 
to statistical fluctuations, given the new number of dofs and 
the number of combinations of states. If this probability was 
low (< 4.6%), then we tried to divide the observations into 
three luminosity states. No source required more than three 
states to be fit acceptably. We list the states in Table [5] For 
each state, we list the observations included in a state (col- 
umn 2), and the combined luminosity and hardness ratios of 
the state (columns 3-6). Unless otherwise noted, we placed 
no requirement on the temporal pattern of the states; LTLVs 
whose states are interspersed temporally may be switching be- 
tween those states or the intrinsic variability may be poorly 
characterized by a simple transition between states. There are 
several categories of luminosity state transitions that indicate 
potentially interesting transitions and a few sources that bear 
individual scrutiny. 

Sources 25, 78, and 110 are SSs that exhibit possible LTLV. 
As SSs, they are more affected by QE degradation than the 
typical source. Among these sources, the first two have a lu- 
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minosity that drops precipitously in the last four observations 
compared to the first observation. Since the last four observa- 
tions all have a similar amount of QE degradation, we reran 
the LTLV detection method only considering the last four ob- 
servations; neither of the sources remained LTLVs. There- 
fore, we do not consider Sources 25 and 78 as LTLVs. Since 
Source 110 does not turn on until the last observation, the QE 
degradation does not complicate its identification as an LTLV. 

Sources 68 and 117 were the only sources broken into three 
luminosity states. The individual luminosities of Source 68 
are qualitatively described by an increase over the first three 
observations, followed by the source turning off, and perhaps 
beginning to turn on again in the last observation. Its two- 
state mode is consistent with this qualitative picture, contain- 
ing a higher-luminosity state of the first three observations 
and a state where the source is not significantly detected in 
the last two observations; however, the probability the two- 
state description is correct is 6. 1 x 10~ 3 . In its best three-state 
mode, Observations 0784 and 4730 are grouped together in a 
lower-luminosity state that is still significantly detected; we 
note that the three-state set where the first two observations 
share a lower-luminosity state, the third observation is the 
higher-luminosity state, and the last two observations are a 
state where the source is not significantly detected has a \ 2 
only 0.46 higher than the best three-state set. The individual 
luminosities of Source 117 vary significantly with no clear 
temporal pattern. In its two-state set, Observation 4729 is in 
the lower-luminosity state and the rest of the observations are 
in the higher-luminosity state; however, the probability the 
two-state description is correct is 4.2 x 10~ 5 . 

Some sources had a lower-luminosity state that was consis- 
tent with no emission at the 3a level, and a higher-luminosity 
state, whose observations were sequential and had emission 
detected at the > 3 a level. We labeled sources following 
this pattern of states as transient candidates. For Sources 68 
and 117, we considered their two-state properties for discus- 
sions of transience. Eleven of the Analysis Sample sources 
(8.9+2g%) are transient candidates (Sources 11, 37, 41, 45, 
51,67,68, 77, 94, 110, and 121). For the eleven transient can- 
didate sources we also calculated the minimum ratio of fluxes 
at the 3a level. Only the flux ratio of Source 110 is above 
10. We believe Source 110 is a truly transient source and in 
§ !4.4l we discuss it in detail. The other transient candidates are 
either transient sources or LTLVs that have a lower luminos- 
ity state that is below the detection limit. All eleven transient 
candidates appeared to either turn on or turn off during the 
4.6 years over which we observed; no source turned on and 
then off, or off and then came back on again. If we assume 
all the LMXBs in NGC 4697 are long-term transients, and 
that our detected transient sources are just those which turned 
off or on during our observation period, then the mean out- 
burst duration for LMXBs in NGC 4697 is 103^ yr. Under 
this hypothesis, the duty cycle of LMXBs would also depend 
on their recurrence timescale. Since our observations find no 
recurring transients, we have adopted the crustal heating re- 
currence tim escale for the Galactic NS-LMXB KS 1731-260 
of ~ 1500 yr (Rutled ge et al.l2002l) . (The adoption of a crustal 
reheating timescale also implies that we are assuming all the 
LMXBs in NGC 4697 are NS-LMXBs. While this is not nec- 
essarily true, we note that the recurrence timescale of Galactic 
long-term transient BH-LMXBs is unknown. As such, we 
present this calculation for illustrative purposes only.) We 
note that this implies a duty cycle of ~ 7%, which is compara- 
ble to the duty cycle obtained comparing active and quiescent 



LMXBs in Galactic GCs (~ 12%; iHeinke etaill200l . We 
also note that the number of transient candidates and possible 
transient candidates in N GC 3379 yields similar duty cycles 
(Brassingt on et al.l2008l) . although their study uses a different 
definition of transience. 

Six of the LTLVs have a > 2a difference in at least one of 
their hardness ratios between their states (Sources 1 1, 45, 1 17, 
121, 134, and 156; Table|5] columns 4-6). These sources ap- 
pear to exhibit luminosity/spectral state transitions. Since the 
lower-luminosity state in Sources 45 and 121 is not even a la 
significant detection, we exclude these sources from this dis- 
cussion. Although Source 1 1 is a Field-LMXB that is a tran- 
sient candidate, it may also be a luminosity/spectral state tran- 
sition source. The luminosity of the lower-luminosity state 
(Observation 4730) is more than 2a significant and appears 
to be 3.4cr softer in H21, than the higher-luminosity state. Al- 
though this softer/fainter to harder/brighter transition is oppo- 
site of the conventional relationship in BH state transitions, it 
is possible that this discrepancy results from the softer bands 
used to calculate the hardness with Chandra (Juett et al. 2008 
in preparation). Alternatively, this source may be transition- 
ing from the thermal state (formerly called high/soft) to the 
steep power law state (formerly called very high). For ex- 
ample, in its 1996-1997 outburst, the Galactic LMXB GRO 
J1655-40 is softer and faint er in its thermal state compare d 
to its steep power law state (Remillard & McClintock 2006). 
Source 1 17, a known AGN, appears to display the more typ- 
ical harder/fainter to softer/brighter transition (e. g., the spec- 
tral c orrelation of Seyfert galaxy NGC 4151 in iPerola et al.1 
119861) : Observation 4729 is ~ 4-5 a harder in i/31 and H i2 
than the two higher luminosity states. Its spectrum seems 
to change most in the 2-6 keV band. Source 134, a sus- 
pected AGN based on its hardness ratios, appears to undergo 
a change in its 1-2 keV band. The lower-luminosity state 
is 2.2a harder in H21, but 2. Oct softer in 7/32 than its higher- 
luminosity state. Finally, source 156, a known foreground 
star, undergoes a transition similar to source 1 1 . The lower- 
luminosity state of source 156 is 2.4cr softer in i/21 than its 
higher-luminosity state. This behavior is not surprising given 
that stellar flares generally heat their corona (Gudel 2004) and 
the higher-luminosity state occurs in the observation where 
this source flares. 

4.3. Long-Term Hardness Ratio Variability 

Our method for identifying LTLV can also be applied to 
hardness ratios; we define sources whose hardness ratios vary 
as Long Term Hardness Variables (LTHVs). Sources that have 
variations in the hardnesses ratios that are significant at the 
2a level have the probabilities that their variations are due to 
statistical fluctuations listed in Table [4] columns 13-15. We 
identify fewer LTHVs than LTLVs because the errors in hard- 
ness ratios are larger than in luminosities. Therefore, we only 
summarize the results for the 83 sources in multiple observa- 
tions and with luminosities detected at the 5a level. We detect 
6.0^|%, 7.3±y%, and 6. sources that exhibit H 2 u #31, 

and H32 variability at the > 2a significance level, respectively. 
Corrections for falsely identified variable sources reduce the 
percentages of hardness variable sources to < 5.3%, < 6.9%, 
and < 5.5%. That is, the numbers of sources selected purely 
based on hardness ratio variations could result from statistical 
fluctuations. 

We also a pplie d the same grouping of observations into 
states as in § 14.21 using individual hardness ratios as opposed 
to luminosities (Table [5] lower portions). Here we discuss 
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sources that were identified as LTHVs in at least two hard- 
ness ratios, or identified as LTHVs in at least one hardness 
ratio and also as LTLVs. Sources 10, 66, and 102 are in the 
first category, while Sources 11, 117, 134, and 156 are in the 
the second. We discuss attempts to resolve the difference in 
their states, and the implications of those attempts. 

For Source 10, the selections by //31 and H32 agree, point- 
ing to a probable decrease in the 2-6 keV count fraction for 
Observation 4729. Although the best-fit states disagree for 
Source 66, matching the H32 selected state to the H31 selected 
state leads to a \ 2 that is only 0.7 higher than the best H32 
selected state. The fraction of 2-6 keV counts is higher for 
Observation 4729. Similarly, matching the H31 selected state 
to the H32 selected state leads to a \ 2 that is only 0.2 higher 
than the best H31 selected state for Source 102; the 2-6 keV 
count fraction is higher for Observations 4727 and 4729. 

Sources 11, 117, 134, and 156 were all identified as lu- 
minosity/spectral state transitions using the LTLV selection. 
The hardness ratio two-state sets in Source 1 17 are consistent 
with the two-state luminosity set; the 2-6 keV count fraction 
is higher in the lower-luminosity state (Observation 4729). 
The H21 and luminosity states match for Source 156, where 
the lower-luminosity state is softer in H21 than the higher- 
luminosity state. Sources 11 and 134 are more difficult to 
resolve as clear luminosity /spectral state transitions as the the 
grouping of observations into states by luminosity and hard- 
ness ratios do not match. If we match the H21 selected state 
to the luminosity selected state for Source 11, the x 2 is 4.1 
higher; however, both sets of states are consistent with the 
lower-luminosity state being softer in H21. The resolution for 
Source 134 is even more problematic. This source appears 
best described by a non-monotonic evolution of hardness with 
luminosity. As the source luminosity increases, it first be- 
comes softer and then becomes harder in all three hardness 
ratio bands; however, the softest observation changes between 
4728 or 4729 depending on the hardness ratio. 

4.4. Source 110: A Transient Supersoft Source 

We consider Source 110 to be a clear transient source. In 
the sum of the first four observations, there were 5 counts in 
the source aperture and 19 counts in the (three times) larger 
background aperture. After correcting for QE degradation 
of a typical source in NGC 4697, the corresponding rate is 
(-1.3 ±3.8) x 10~ 5 cnts _1 ; the source is clearly undetected. 
We also calculated the corrected rate including observations 
4727-4729, which have a similar level of QE degradation, to 
be (-1.0±4.6) x 10~ 5 cnts _1 . In the final observation, there 
were 87 counts in the source aperture and 4 counts in the back- 
ground aperture, corresponding to a QE degradation corrected 
rate of (32.6 ±4.0) x 10~ 4 cnts _1 . The ratio of the count rate 
in Observation 4730 to the other observations is > 27 at the 
3cr level. Therefore, we only consider data from Observation 
4730 hereafter. Within the 0.3-6.0keV band, only two of the 
source aperture photons were above 1 keV (1.03 and 3.9keV). 
Source 110 is clearly supersoft. Its 3a limits on observed 
hardness ratios are //21 < -0.89 and H31 < -0.88. Analysis 
of our flanking-field HST-ACS observations indicates Source 
110 is not associated with an optical source, although it is 
about 1"9 away from a GC candidate. From simulated point 
sources at the location of Source 110 we placed 90% confi- 
dence limits on magnitudes for a single HST-ACS image of 
g > 26.3 and z > 25.3 in AB mag. Since Source 110 is in 
an overlapping region of two fields, these magnitudes corre- 
spond to the 99% confidence limit on an optical counterpart 



to Source 110. 

Since Source 1 10 is both transient and SS, we explored its 
spectrum in greater detail. We examined both the raw spec- 
trum and a binned spectrum. Since Source 110 has only 87 
gross counts in Observation 4730, we created five bins with 
at least 16 counts in the 0.3-1 keV range and a bin with two 
counts from l-6keV. These bins were weighted by Gehrels 
errors. 

Luminous SSs in our Galaxy and M31 are generally be- 
lieved to be accreting WDs. Such sources are often fit by 
a blackbody (bb) model. Compared to the best-fit for this 
model, the spectrum appears to rise too quickly within the 
0.3-1 keV band to be consistent with Galactic absorption in 
the raw spectrum. In the binned spectrum, the best-fit model 
with Galactic absorption was strongly rejected, x 2 = 16.78 for 
4 dof (kT bb = 129 ± 16eV). This mo del is also hotter than typ- 
ical fi ts to SS-WDs (15-80eV; iKahabka & van den Heuvel 
1997). If we allow the absorption to vary, we do get a 
high absorbing column and lower temperature, Nh = 1 -2*q g x 
10 22 cm 2 and kT bb = 55^ eV (90% individual confidence lim- 
its), for an acceptable fit with \ 2 = 2 04 for 3 dof. How- 
ever, the implied X-ray and bolometric luminosity in the 0.3- 
lOkeV, L x = 3.0 x lO^ergss" 1 and L bol = 1.6 x lO^ergss" 1 
are much too large for a WD at the distance of NGC 4697. 
Even the lowest X-ray and bolometric luminosities allowed 
by the uncertainties at 90% confidence for jointly varying ab- 
sorption and temperature, L% = 1.6 X 10 39 ergss _1 and L bo \ = 
3.0 x 10 39 ergss _1 , are much too luminous for a SS-WD 
in NGC 4697. These luminosities come at hotter tempera- 
tures, 90eV, than seen in SS-WDs. Although more detailed 
WD atmosphere models can lead to luminosities lower by 
an order of magnitude than derived from the bb mode l (e.g., 
iHeise et al.|[l994T: iBalman et al.lll998t lOrio et al.ll2007l) . they 
also tend to increase the temperature. Thus, we consider it un- 
likely that Source 1 10 is a SS-WD in NGC 4697. On the other 
hand, our optical limit suggests that it is unlikely that Source 
110 is a SS-WD within the Milky Way. Of the likely com- 
panions for a SS-WD, main sequence stars earlier than F5, 
red giants overflowing the ir Roche lobe, and asymptotic gian t 
branch stars with winds (Kahabka & van den Heuvellll997l) . 
an F 5 main-sequence star is the faintest with My ^3.5 (pg. 
107 iBinney & Merrifieldl fl987h . We would detect such a 
source out to ~ 300kpc in our HST images. Thus, Source 
1 10 is unlikely to be a SS-WD, either as a foreground Galac- 
tic source or a source in NGC 4697. 

A very soft, luminous X-ray transient with an excess ab- 
so rption column has been o bserved in M101 (M101 ULX- 
1: iKong & Di Stefanoll2005h . The combination of spectra, 
outburst luminosity, and transience was used to suggest that 
source is an intermediate-mas s BH (IMBH). The disk black- 
body model (hereafter diskbb, Mitsuda et al. 1984) is a multi- 
color blackbody model that is often used to represent BH 
soft X-ray transients (SXT s), which have kT^^ < l-2keV 
dTanaka & Shibazaki 1 19961) . As with the simple blackbody 
model, the spectrum is poorly fit by a model with Galac- 
tic absorption; kT^^ = 163 ± 27 eV yields % 2 = 20.05 for 4 
dof. If we allow the absorption to vary, the absorption and 
temperature are strongly correlated and span a fairly large 
range, but the fit is much improved with a \ 2 = 2.02 for 3 
dof. We find N H = 1.2^ x 10 22 cm 2 and kT diskbb = 59!geV 
(90% individual confidence limits). The 0.3-10keV X-ray 
luminosity of acceptable spectral models (within the 90% 
two-dimensional confidence interval for varying absorption 
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and temperature) ranges from 3.5 x 10 - 6.8 x 10 ergs s . 
The bolometric luminosity range is even larger, 1.4 x 10 40 < 
^boi < 4.8 x 10 48 ergss _1 . We use the normalization of the 
di skbb model to estim ate a BH ma ss following equat ion (3) 
in IFabian et all (120041) . Following iFabian et all (120041) . we 
also assume a correction factor of /o = 1.35, an inclination 
angle of i = 60°, and an innermost disk radius of 6GMbh/c 2 , 
where Mbh is the BH mass. Dividing our derived bolomet- 
ric luminosity by the Eddington luminosity, Z^dd = 1-26 x 
10 38 (MBH/A^o)ergss _1 , gives the Eddington efficiency of the 
accretion disk, which we use later in evaluating the plausibil- 
ity of the model. The value of the disk blackbody temperature 
also provides a second estimate of the black hole mass, which 
is rou ghly M Bii /M G ~ (kT diskbb / 1 ke V)" 4 (Makishim a etalj 
2000). Considering models within the 90% two-dimensional 
confidence interval for varying absorption and temperature, 
these two estimates of BH mass agree to within a factor of 
two for absorption columns between 0.5 and 1 x 10 22 cm 2 . 
Within this region, our source is consistent with BH masses 
of - (7.2 x 10 3 — 1.1 x 10 6 )M Q accreting at - 1.6-4.7% of the 
Eddington limit. Thus, this source is consistent with a rather 
massive IMBH or a small supermassive black hole about a 
factor of 2-3 s maller in mass th an the one at the center of 
the Milky Way (Ghe z et alj|2000h . Understanding the origin 
of such a large BH outside of the center of a galaxy or of a 
GC is a formidable challenge to our current theories of BH 
formation. 

Some AGNs have ul trasoft spectra. For example, 
iPuchnarewicz et alJ (Il992h identified 53 ultrasoft AGN can- 
didates with Einstein Observatory. These sources had 0.16- 
0.56keV count rates nearly three times that of their 0.56- 
1.08keV count rates. This corresponds to steep power-law 
photon indices, T > 3, with Galactic abso rption. In a s t udy of 
variability in soft X-ray selected AGNs, iGrupe et ail (2001) 
found that AGNs with steeper X-ray spectra had larger short- 
term variability than those with flatter spectra. Two of the 
transient sources, IC 3599 and HE 0036-5133 (also called 
WPVS 007) were also previously reported as ultrasoft. These 
sources may be analogs to Source 110. A variety of compli- 
cated models have been used to explain the spectra of AGNs 
with excess soft emission compar ed to the power-law mod - 
els derived from harder bands (e.g., Gierlihsk i & Dond2004l) . 
Given the relative paucity of our data on Source 1 10, we chose 
to only consider simple power-law and blackbody models. 
The power-law model with Galactic absorption, T = 4. 23^0 ^, 
is an exceptionally poor fit, \ 2 = 30.20 for 4 dof. Even with 
a free absorption, power-law models with reasonable photon 
indices, T < 10, did not fit our data well, \ 2 > 12.25 for 3 
dof. From above, we know that the disk blackbody model 
is a good statistical representation of the spectra. However, 
our derivation of the BH masses in the above paragraph was 
predicated on the source being at the distance of NGC 4697. 
As we allow the source to lie farther behind NGC 4697, the 
agreement between the mass derived from the normalization 
and the mass derived from the temperature begins to break 
down. Thus, it is unlikely that the disk blackbody model 
fully represents the physical mechanism behind the emission. 
Regardless of the model we choose to represent the emis- 
sion, we measure an absorbed X-ray flux (0.5-8.0keV) of 
Fx ~ 7 X 10~ 15 ergss _1 cm" 2 . At that flux level, we only ex- 
pect — 3.6 sources unrelated to NGC 4697 for galactocen- 
tric radii interior to that of Source 110. Thus, the probability 
Source 1 10 is a background AGN is — 3.3% from the X-ray 



data alone. If we assume that the transient behavior observed 
was not an extreme flare that began between 2004 February 
02 and August 18, but ended before the HST-ACS flanking 
field observations began in 2005 December 19, we can place 
an additional constraint using the logarithm of its X-ray-to- 
optical flux ratio, log(,Fx /F opt ), where F opt is the optical flux 
in the z band. Based on our non-detection of a counterpart 
to Source 1 10, its log(F x /F opt ) is greater than 0.9 at the 99% 
confidence level. Such a large ratio is expected only from 
the most obscured AGNs. There are two problems with sug- 
gesting Source 110 is an obscured AGN. First, the most ob- 
scured AGNs can have columns one to two orders of magni- 
tude larger than the estimated column for Source 110. In fact, 
its column is actually near the division point ty pically used 
for obscured and unobscured AGN (10 22 cm 2 : iTreister et alJ 
2004). Second, very few sources are predicted to have su ch 
high X-ray-to-optical flux ratios. In Treist eret alJ ([2004), a 
model based on the AGN unification paradigm predicts only 
~ 13% of AGNs will have log(F x /F opt ) > 0.8. In summary, 
although some of its X-ray properties may be consistent with 
transient ultrasoft AGN behavior, the combination of the ex- 
pected number of background AGNs with the high X-ray-to- 
optical flux ratios makes it unlikely (probability — 0.4%) that 
Source 1 10 is a transient ultrasoft AGN. 

5. CONCLUSIONS 

Our multi-epoch Chandra observations were planned to 
shed light on LMXB variability in NGC 4697, the nearest, 
optically luminous, elliptical galaxy. We measure variability 
on timescales of seconds to years. Although we do not de- 
tect any periodic LMXBs at significant levels, we detect five 
sources with significant variability on timescales of seconds 
to hours (Source 57, 71, 75, 155, and 156). 

When we compare the luminosities of sources between ev- 
ery pair of observations, we determine that 6.7 ± 1.0% of 
the Analysis Sample sources are variable over timescales of 
— 10-2000 days. There is a slight increase in the frac- 
tion of variable sources with the variability timescale, but 
it is not found to be statistically significant. Our number 
is roughly consistent with that fo und for Centaur us A for 
a variability timescale of 164d dKraft et alJ 1200 lh and for 
NGC 4636 for a day timescale, but i s lower than that found fo r 
a 1 1 15 d timescale in NGC 4636 dPosson-Brown et alj|2008l) . 
New multi-epoch observations of NGC 1023, NGC 3379, 
NGC 4278, NGC 4365, and Centaurus A will help fill in 
more variability timescales. As this phase space is popu- 
lated, we will be able to test if more sources are variable on 
longer timescales than shorter timescales. Increased knowl- 
edge of the temporal evolution of variability will impact mod- 
els that explore the nature of the observed LMXBs (e.g., 
iPiro & Bildstenll2002l) . 

We are more sensitive to variability when we consider all 
five observations simultaneously. We find that w 16 ±4% of 
Analysis Sample sources do not have a constant luminosity 
over all five observations. If we apply a si milar requirement 
for a s ource to be considered variable to the Brassingto n et alJ 
(2008) analysis of LMXBs in NGC 3379 as in this study, we 
find that the number of variable sources is 12 ±5%, consis- 
tent with what we find for NGC 4697. The fraction of GC- 
LMXBs that are variable is statistically indistinguishable from 
the fraction of variable Field-LMXBs in our sample due to the 
small number of variable sources. To better test for any dif- 
ference between the two populations, multi-epoch variability 
from multiple galaxies will need to be combined. 
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For the detected sources we group observations into sets 
of luminosity and hardness ratio states. Among these variable 
sources, we identify eleven transient candidates. If we assume 
all the LMXBs in NGC 4697 are long-term transients, we es- 
timate that they are on for ~ lOOyr and have a 7% duty cycle. 
The latter is comparable to the duty cycle obtain ed comparing 
active and quiescent LMXBs in Galactic GCs (Hein keet alj 
l2003h . 

The mean outburst duration for LMXBs in NGC 4697 of 
~ lOOyr is consistent with a lower limit of ~ 50 yr found 
for 33 luminous LMXBs in NGC 1399 and M87 that ha d 
X-ray luminosities exceeding 8 x 10 38 ergs s" 1 (Ilrwinll2006h . 
These 33 sources were assumed to be BH-LMXBs because of 
their sustained super-Eddington luminosities over several year 
timescales, and most likely represent extragalactic analogs of 
the Milky Way X-ray binary GRS1915+105. Since a major- 
ity of the sources in NGC 4697 are probably NS-LMXBs, 
there does not appear to be a measurable difference between 
the outburst durations of long-term transient NS-LMXBs and 
BH-LMXBs. 

We also find four sources with a combined luminos- 
ity/spectral transition, only one of which is likely associated 
with NGC 4697. One limit to our ability to detect such 
sources is the duration of our single observations. Although 
observations of NGC 3379 and NGC 4278 have longer single 
observation durations, this issue is likely to be best addressed 
by a recently completed set of six, 100 ks observations of Cen- 



taurus A, which is at least 2.5 times closer than any other mas- 
sive early-type galaxy. 

Source 110 is clearly a supersoft transient source with no 
detected counterpart in HST observations. Although it only 
has « 85 X-ray counts, our spectral analysis rules out the pos- 
sibility that this source is a luminous supersoft WD. The high 
X-ray-to-optical flux ratio and relative small number of ex- 
pected background AGNs argue against the source being a 
transient ultrasoft AGN. Our spectral fit could be consistent 
with an accreting BH with a mass in the range of ~(7 x 10 3 - 
10 6 )M Q ; however, the origin of such a large BH outside of the 
galaxy center or a GC is a formidable challenge to our current 
theories of BH formation. Further observations of NGC 4697 
could increase the information we can extract from the spec- 
tra of Source 110; however, more work must also be done to 
increase our understanding of the X-ray spectra of non-stellar- 
mass black holes. 
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States of Variable X-ray Sources 
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